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Abstract. Under certain plausible assumptions, M. Rubinstein and P. Sarnak solved 
the Shanks-Renyi race problem, by showing that the set of real numbers x >2 such 
that T:{x;q,ai) > iT{x;q,a2) > ■ •• > ^{x'jq^ar) has a positive logarithmic density 
^q;ai,...,ar- Furthermore, they established that if r is fixed, Sq-ai,...,ar — ^ 1/r! as g oo. 
In this paper, we investigate the size of these densities when the number of contes- 
tants r tends to infinity with q. In particular, we deduce a strong form of a recent 
conjecture of A. Feuerverger and G. Martin which states that Sq-ai,....ar. — o(l) in this 
case. Among our results, we prove that Sq-ai,...,ar ~ l/*"' in the region r — o{^/\ogq) 
as (7 — )■ oo. We also bound the order of magnitude of these densities beyond this range 
of r. For example, we show that when logg < r < (j>{q), Sq-ai,...,ar. -^e g"^"*"'. 



1. Introduction 

A classical problem in analytic number theory is the so-called "Shanks-Renyi prime 
number race" which concerns the distribution of prime numbers in arithmetic progres- 
sions. As colorfully described by Knapowski and Turan in [TT], let g > 3 and 2 < r < 
0(g) be positive integers, and denote by Ar{q) the set of ordered r-tuples of distinct 
residue classes (oi, 02, ... , a,.) modulo q which are coprime to q. For (oi, 02, ... , a^) G 
Ar{q), consider a game with r players called "1" through "r", where at time x, the 
player "j" has a score of 7c{x;q,aj) (where 7r(a;;g,a) denotes the number of primes 
p < X with p = a mod q). As x — )■ 00, will all r! orderings of the players occur for 
infinitely many integers x? 

It is generally believed that the answer to this question is yes for all q and all 
(ai, 02, . . . , Or) G Ar{q). An old result of Littlewood ^J] shows that this is indeed 
true in the special cases (g, 01,02) = (4,1,3) and (g, 01,02) = (3,1,2). Since then, 
this problem has been extensively studied by many authors, including Knapowski and 
Turan [TT] , Bays and Hudson [T] and [2] , Kaczorowski |S] , P and [TU] , Feuerverger and 
Martin [1] , Martin [15] , Ford and Konyagin [6] and [7] , Fiorilli and Martin [5] , and the 
author [121 and 
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A major breakthrough was made in 1994 by Rubinstein and Sarnak who completely 
solved this problem in [16], conditionally on the two following assumptions: 

• The Generalized Riemann Hypothesis (GRH): all nontrivial zeros of Dirichlet 
L-functions have real part equal 1/2. 

• The Linear Independence Hypothesis (LI) (also known as the Grand Simplicity 
Hypothesis): the nonnegative imaginary parts of the nontrivial zeros of Dirichlet 
L-functions attached to primitive characters are linearly independent over Q. 

Rubinstein and Sarnak proved, under these two hypotheses, the stronger result that 
for any (ai, . . . , a^) G Ar{q), the set of real numbers x >2 such that 

n{x; q, ai) > tc{x; g, 02) > ■ ■ ■ > tt{x; q, a^), 

has a positive logarithmic density, which shall be denoted throughout this paper by 
Sq;ai,...,ar- (Recall that the logarithmic density of a subset S of M is defined as 

6s := lim [ ^, 

x^oo logx Jt^sn[2,x] t 

provided that this limit exists). To establish this result, they constructed an absolutely 
continuous measure fiq-ai,...,ar for which 



(!•!) 6q-ai^,,,^ar / '^/^q;ai,...,ar ("^l; • • • ; -^r) ■ 

J Xl>X2>--->Xr 

Among the results they derived on these densities, Rubinstein and Sarnak showed 
that in an r-way race with r fixed, all biases disappear when g — )■ 00. More specifically 
they proved 

(1-2) lim max |r!(5g;ai,. 1| = 0. 

q-S>oo {ai,...,ar)&Ar{q) 

Recently, Fiorilli and Martin [5] established an asymptotic formula for the density in a 
two-way race, which allows them to determine the exact rate at which 6q-ai,a2 converges 
to 1/2 as q grows. Shortly after, the author [il2j succeeded to obtain an asymptotic 
formula for Sq.ai,...,ar foi' ^"^J fixed r > 3 as g — )■ cxd, in which the rate of convergence to 
1/r! is surprisingly different from the case r = 2. 

However, as far as the author of the present paper knows, no results have been 
obtained on the size of the densities Sq-ai,...,ar if r — )■ 00 as g — )■ 00. In Feuerverger 
and Martin conjectured that in this case we should have Sq-ai,...,ar = o{l). They also 
asked whether one can prove a uniform version of the result of Rubinstein and Sarnak 
(1.2), namely that this statement holds in a certain range r < ro(g) for some ro(g) — > 00 
as g — 7- 00. 

Conjecture 1.1 (Feuerverger-Martin) . We have 

lim max So■a^ a = 0, 

g-s-oo {ai,...,ar)eAr(q) 
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for any arbitrary function r = r{q) tending to infinity with q. 

In the present paper, we investigate the order of magnitude of 5Q;ai,...,a^ when the 
number of contestants r — )■ oo as g — > oo. In particular, answering the question of 
Feuerverger and Martin, we establish a uniform version of (1.2), and obtain a strong 
quantitative form of Conjecture 1.1. 

Theorem 1.1. Assvme GRH and LI. Let q he a large positive integer. Then, for any 
integer r such that 2 < r < \/\og q we have 

If f r'^ 



(1-3) 3q;ai,-,ar ~ ^> 



r! \ \logQ' 
uniformly for all r-tuples (oi, . . . , a^) G Ar{q)- 

As a consequence. Theorem 1.1 implies that (1.2) holds true in the range r = 
o{y/\ogq) as g — )■ oo. Indeed in this region of r, all biases disappear when g — )■ oo, 
namely 

1 

H 

uniformly for all r-tuples (oi, . . . , a^) G Ar{q)- Moreover, one can also deduce that if 
Co > is a suitably small constant and r < co^/\ogq, then uniformly for all r-tuples 
(ai, . . . , a^) e Ar{q) we have 

(1.4) Sq;ai,...,ar ^ ^• 

Note that 1/r! = exp(— r logr + r + O(logr)) by Stirling's formula. Our next result 
shows that the densities Sq-ai,...,ar have roughly the same asymptotic decay in the range 
\/log q <C r < (1 — e) log q/ log log q, for any e > 0. 

Theorem 1.2. Assume GRH and LL For any e > 0, if q is large and ^/\ogq <^ r < 
(1 — e) logg/ log log g is an integer, then 

Sqm,-,ar = exp l^-rlogr + r + O ( logr + 

uniformly for all r-tuples (ai, . . . , a^) € Ar{q)- 



^2 



logg 



It would be interesting to determine the order of magnitude of the densities Sg.ai,...,ar 
beyond the region r < (1 — e) logg/ log log g. Unfortunately, this range seems to be the 
limit of what can be achieved using our method. Nevertheless, we can use Theorem 1.2 
to obtain an upper bound for Sg.ai,...,ar beyond this range of r. 

Theorem 1.3. Assume GRH and LL For any e > 0, ifq is large and (1— e/2) log q/ log log 

r < (f>{q) is an integer, then 

1 

max dg.ai,...,ar <e -r-. 

{ai,...,ar)eAr{q) 
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The paper is organized as follows. In Section 2, following the work of Rubinstein 
and Sarnak, we shall construct the measure fiq;ai,...,ar as a probability distribution 
corresponding to a certain random vector and study its covariance matrix and large 
deviations. In Section 3, we investigate the Fourier transform of fiq-ai,...,ar and show 
that in a certain range fiq;ai,...,ar can be approximated by the Fourier transform of a 
multivariate normal distribution having the same covariance matrix. In Section 4, we 
study properties of multivariate normal distributions and prove Theorems 1.1, 1.2 and 
1.3. 

2. The measure Hq-ai,...,ar 

We begin by developing the necessary notation to construct the measure fiq-ai,...,ar, 
following the work of Rubinstein and Sarnak [16j . For (oi, a2, . . . , Or ) G Ar{q) we in- 
troduce the vector-valued function 

Eq-ai,...,ar{x) := {E{x; q, Gi) , . . . , E{x; q, Qr)) , 

where 

E{x; g, a) := (0(g)7r(x; q, a) - 7r(x)) . 

The normalization is such that, if we assume GRH, Eq.ai,...,ar{^) varies roughly bound- 
edly as x varies. Moreover, for a nontrivial character x modulo q, we denote by {7^} the 
sequence of imaginary parts of the nontrivial zeros of L{s, x)- Let Xo denote the princi- 
pal character modulo q and define S = ^x¥=xo ^od q{lx} ■ Furthermore, let {^(7x)}7x6S' 
be a sequence of independent random variables uniformly distributed on the unit circle. 

Rubinstein and Sarnak established, under GRH and LI, that the vector-valued func- 
tion Eq.ai,...,ar Las a limiting distribution /ig;ai,...,a^, where f^g-ai,...,ar is the probability 
measure corresponding to the random vector 

Xq-a^,...,ar = (^(?, ^l), . . . , X {q, O^)), 

where 

X(,.a) = -C,(a)+ ^2Re(x(a)[/fa))_ 

Xt^Xo 7x>0 \/i +7x 

X mod q ^ 

and 

Cq{a) :=-!+ Yl 1- 

b^=a mod q 
l<b<q 

Note that for (a, q) = 1 the function Cq{a) takes only two values: Cq{a) = —1 if a is a 
non-square modulo q, and Cq{a) = Cq{l) if a is a square modulo q. Furthermore, an 
elementary argument shows that Cq{a) < d{q) q^ for any e > 0, where d{q) = 1 
is the usual divisor function. 
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To investigate the distribution of the random vector Xq.^ai,...,ar we shall first com- 
pute its covariance matrix Covq-ai,...,ar (the covariance matrix generalizes the notion of 
variance to multiple dimensions). Recall that the j,k entry of the covariance matrix 
corresponds to the covariance between the j-th and A;-th entry of the random vector. 



Lemma 2.1. The entries of CoVq.ai,...,ar 

C0Yq.au...,ar{3,k) 



fVar(g) if j = k 

]^Bg{aj,ak) ifjj^k, 



where 



VarW:=2 $3 ^ and B,(a, 6) J^^M^^- 

XJ^XO 7x>0 4 'X Xt^XO 7x>0 ^ 

X mod q X mod q 

Proof First, note that E(X(g, a)) = — Cg(a) since E(t/(7x)) ~ ^ 7x- Therefore, 

Gov,, 

;ai,...,ar (j? ^) Gquals 

E((X(g, a,) + C,(a,)) (X(g, a^) + C,(afc)) 



E 



/ ^ \ 

2-^ 2^ 2^ 2^ /i 7 /i — ^ 

. Xt^XO 7x>0 ^7^X0 7v.>0 V4+7x\/4+7v. , 

\ X mod (jf i/" mod q V V / 



Since E(?7(7x)f/(7^)) = for all 7^5 7i/; and 



1 if X = and 7x = 7^, 
otherwise. 



we deduce that 

Xt^XO 7x>0 4 + 7x 

X mod Q 

which implies the result. □ 

Our next lemma gives the asymptotic behavior of Var(g) along with the maximal 
order of Bq{aj, a^). This was established in [12], and we should also note that it follows 
implicitly from the results of [5]. 

Lemma 2.2. Assume GRH. Then 

(2.1) Var(g) = 0(g) log g + O(0(g) log log g), 
and 

(2.2) max i?,(a, 6) X 0(g). 

(a,fe)e^2(g) 
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Proof. First, the asymptotic formula (2.1) is proved in Lemma 3.1 of [12]. Now, the 
fact that Bq{aj,ak) ^ (t>{q) is proved in Corollary 5.4 of [12j, while Proposition 5.1 of 
implies Bg{a, —a) ^ (t>{q)- □ 



Here and throughout we shall use the notations ||t|| = ■\^Yl]=i^'j ^^"^ l^loo = 
maxi<j<r \tj\ for the Euclidean norm and the maximum norm of t G M*" respectively. 
Our next result is an upper bound for the tail of the distribution fiq-ai,...,ar- This was 
established in Proposition 4.1 of [12] in the case where r is fixed. 

Lemma 2.3. Let q be large and 2 < r < 0(g) be a positive integer. Then for R > 
a/ (f){q) logg we have 

/x,;a......(|xU>i?)<2rexp (-4^(|l^), 

uniformly for all (ai, . . . , a.,.) G Ar{q). 
Proof. First, we have 

r r 

/Xg;ai,...,a.(|x|oo > i?) = P(|X,;,,,...,aJoo > R) < J] P(X(g, a,) > /?)+^ P(X(g, a,) < -R). 

i=i i=i 
We shall bound only P(X(g, aj) > R), since the corresponding bound for P(X(g, aj) < 
—R) can be obtained similarly. Let s > and (a, g) = 1. Then we have 



X mod <? 



Xt^XO 7x>0 
X mod (jf 




where Io{t) '■= 'Yl^=o{^ / '^)'^^ / ^^-^ ^'^^ modified Bessel function of order 0. Hence, using 
the Chernoff bound along with the fact that Io{s) < exp(s^/4) for all s G M we derive 

P(X(g,a) > R)< e-^^E (e^^^^''^)) < exp (^sR - sCq{a) + yVar(g)^ . 

The lemma follows upon choosing s = R/{(f){q) logg), since Cg{a) = q°^^^ and Var(g) ~ 
(j){q) logg by Lemma 2.2. □ 

3. The Fourier transform jlq-ai,...,ar 

Throughout the remaining part of the paper we shall assume both GRH and LI. 
Moreover, we will use the following normalization for the Fourier transform of an inte- 
grable function / : — C 

/(tl, . . . , tn) = / e-*(*^^^+-+*"^")/(a;i, ■ ■ ■ , Xn)dx, . . . dXn. 
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Then if / is integrable on R" we have the Fourier inversion formula 

/(xi, . . . , x„) = (27r)-" [ e^(*i-i+-+*"-")/(t,, . . . , Qdt, . . . dt^. 
Similarly we write 



for the Fourier transform of a finite measure i/ on M". 

Rubinstein and Sarnak [16] established the following explicit formula for the Fourier 
transform of fig 

;ai ,...,ar 



(3.1) fig 

;ai ,...,ar 



X¥=X0 7x>0 
X mod 9 



i +^2 
4 ^ 'X 



where Jo(2;) = Sm=o(~-'-)™('^/^)^'^/"^'^ Bessel function of order 0. 

Our first result shows that in the range ||t|| < Vai{q)^^^'^^°^^\ the Fourier transform 
fiq-ai,...,ar{'tiy ■ ■ ■ jtr) IS Very closc to the Fourier transform of a multivariate normal 
distribution whose covariance matrix equals CoVq-ai,...,ar- 

Proposition 3.1. Let q be large, 2 < r < logg be a positive integer, and (ai, . . . , ar) G 
Ar{q)- Then in the range ||t|| < Var(g)^^/^ log^ g we have 



flg;ai,...,aritl, . . . , tr) = exp ( - - 1 CoVq.ai,...,ar't ) ( 1 + O 



d{q) log^ q 



Proof. First, the explicit formula (3.1) yields 



log flq.ai,...,ar{tu tr) = ^ ^ log Jq 



x¥=xo 7x>0 

X mod q 



4 ^ 'X 



+0 I iitii^ic 



Using Lemma 2.2 along with the standard estimate 0(g) ^ q/ log log g, we deduce that 
the error term above is ^ q~^^'^d{q) log^ g. On the other hand note that 



ELiX(aj)ij 



4 ^ 'X 



< rlltll < 1 



if g is large enough. Hence, using that log Jo{z) = —2^/4 + O(z^) for \z\ < 1 we obtain 



log/ig;ai,...,a,(^l, • • • ,tr) 



(3.2) 



=xa 
od 

/ 



Ej=ix(%)ij 



Xt^XO 7x>0 
X mod q 



4 ^ 'X 



+ 



4||j.||4 

r t 



v 



X^XO 7x>0 ^4 ^ ^X^ 
X mod q 



+ 



d{q) log ^ g 
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Since ^ ^^^^^^ S7^>o V(i + 7x)^ ^ Var(g), it follows that the error term in the above 



X mod q 



estimate is <^ g ^^'^d{q) log^g. On the other hand, the main term on the RHS of (3.2) 
equals 



XJ^XO 7x>0 4 'X l<j,k<r 
X mod q 



l<j,fc<r 



by Lemma 2.1. 



□ 



Next, we show that fiq;ai,...,ar{'ti, ■ ■ ■ ,tr) is rapidly decreasing in the range ||t|| > 
Var(g)~^/^. In particular, the following result is a refinement of Proposition 3.2 of 
which takes into account the dependence of the upper bounds on r. 



Proposition 3.2. There exists a constant ci > such that, if q is large and 2 < r < 
Ci log q, then uniformly for all ( ) G Ar{q) we have 



exp 



0(9) 



tl 



\P^q;a^,...,ar{tl,---,tr)\ < { exp 



exp 



(logg)^ 
(j){q) \ogq. 



if\\t\\ > 400, 
^/(logg)-2 < ||t|| < 400, 
^/||t|| < (logg)-2. 



Before proving this result we first require the following lemma. 

Lemma 3.3. Let q be large and 2 < r < (j){q)/A be an integer. For a 
Ar{q) and t G W we denote by Mg^^i't) the set of nontrivial characters x ^od q such 
that 



EUxMt-i > l|t||/2. Then 



|M„(t)| > 



Proof. Let 



x¥=xo 

X mod q 



E 

X mod q 



(3.3) r r / ^ \^ ^ /r^2 

= YY^^^'' Y x{a,)Wk) - [Y^n = ^ii)Y^' - [Y^^ 

j=l k=l xmod? Vj=l / j=l Vi=l 

>(</>(?) -r-)||t|r, 
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by the Cauchy-Schwarz inequality. Therefore, using that 
rlltlP, we deduce 



E;=ix(a,)t,| <(E;=iit,i) < 



xeAf,,a(t) 



X^A^g,a(t) 



<r|Mg,a(t)|||tf + 



0(g) 



Combining this estimate with (3.3) completes the proof. 



□ 



Proof of Proposition 3.2. First, assume that ||t|| > 400. For any nontrivial character 
X mod q we define 



F{^.x) ■■= n ^0 

7x>0 



2x 



i + ^2 
4 ^ 'A 



Then, it follows from Lemma 2.16 of [5] that 

(3.4) |F(a:,x)F(x,x)| <e-^ 

for X > 200. Moreover, the explicit formula (3.1) implies 



\P'q;ai,...,ar (^1; • • • 5 



n 

x¥=xo 

X mod q 



If X e Mg,a(t) then 



X]y=i x('2j)^j ^ 200. Furthermore, note that x ^ ^g,a(t) if and 



only if X £ ^g,a(t). Hence, using (3.4) along with the trivial bound \F{x, x)| < 1 (since 

|<^o(a^)| < 1) we derive 

(3.5) 



Ag;ai,...,a^(^l, • • • , ^r) P < Y\. 




r 


xeAf,,a(t) 






= n 




r 


xeM,,a(t) 







X 



,X]F 



< exp 



E 

X6M5,a(t) 



<exp(--|M,,.(t)|||t| 



since every character in a(t) appears once as x and once as x in the product on the 
RHS of (3.7). Combining this inequality with Lemma 3.3 yield the desired bound on 

|/ig;ai,...,a,(tl, • • • ,^r)| m this CaSC. 

Let e = (logg)~^ and suppose that e < < 400. If x ^ ^g,a(t) then 



4 ^ 'X 



> 



i + ^2 
4 ^ 'X 
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We also note that if q is sufficiently large then e (| + 7^) < 2e < 1. Therefore, since 
Jo is a positive decreasing function on [0, 1] and |Jo(2;)| < Jo(l) for all z > 1, we get 



|/ig;ai,...,a,(tl,...,tr)| < JJ 11 

X6M,,a{t)7x>0 



Jo 



1+^2 
4 ^ 'x, 



Furthermore, using the standard bound |Jo(a^)| < exp(— a;^/4) for |a;| < 1, we deduce 
that 



(3.6) 



e 1 

|/ig;ai,...,a.(tl,...,t.)| < exp I -- J] 1^1—2 

XeM,,a(t)7x>0 4 + ^X 



Let A^(T, x) denote the number of 7^ in the interval [0, T]. Then, we have the classical 
estimate (see Chapiters 15 and 16 of [3]) 

iV(T,x) = ^ log |^ + 0(loggT), 

where q* is the conductor of x- Hence, if T = log^ q then N{T, x) ^ log^ q- This yields 

2^ 1 + ^2 - 1+^2 log-^ o ■ 

7x>0 4 ^ 'X 0<7x<log2 9 ^ '^^ ^ ^ 

The upper bound on |/ig;ai,...,a,. (^i, • • • ,^r)| then follows upon inserting this estimate in 
(3.6) and using Lemma 3.3. 

Finally assume that ||t|| < (logg)-2. If q is large enough then 



1+^2 
4 ^ 'X 



< rlltll < 1. 



Hence, using that | Jo(3;)| < exp(— a;^/4) for |x| < 1 we obtain from the explicit formula 
(3.1) 



(3.7) 



|/iq;ai,...,a,(tl, • • • , ^r) | < exp 



E E 



\T.]=iX{aj)tj 



, Xt^XO 7x>0 
y X mod q 



4 ^ 'X 



Furthermore, Lemma 2.2 yields 



E E- 

Xt^XO 7x>0 
X mod q 



ELiX(ai)ij 



i + ^2 
4 ^ 'X 



5Z 5Z T— 5Z X(ai)x(afc)^j4 



+ ^2 

Xt^XO 7x>0 4 'X l<j,fc<r 



X mod g 

Var(g) 



(t? + --- + t^)+ Bg{aj,ak)tjtk 

i<j<k<r 
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Since 



J2 \tM< <r\\tr, 

i<j<k<r \j=l / 

by the Cauchy-Schwarz inequality. Thus, if r < ci log q where ci > is suitably small, 
then 



ELix(ai)^j 



2 



- + 



X mod g 

Inserting this estimate in (3.7) completes the proof. □ 

4. The asymptotic behavior of the densities 5g-ai,...,ar'- Proof of 

Theorems 1.1, 1.2 and 1.3 

We showed in the previous section that in a small region around 0, the Fourier 
transform of fiq-ai,...,ar can be approximated by the Fourier transform of a multivariate 
normal distribution whose covariance matrix equals CoVq-ai,...,ar- If wc normalize by 
^yVa.r(q) then Proposition 3.1 above implies that in the range ||t|| < log^ q we have 
(4.1) 

l^q;ai,...,ar , , • • • , , , , = exp --t Ct [1 + 0' 



where C is an r x r symmetric matrix whose entries arc 

\ Var(g) logg 

Let A^r (e) denote the set of r x r symmetric matrices A = (ajk) such that ajj — 1 for 
all 1 < i < r and \ajk\ < e for all 1 < j 7^ /c < r. In order to prove Theorems 1.1-1.3, we 
need to investigate multivariate normal distributions whose covariance matrices belong 
to M.r{€) where e <^ 1/logg' is small. To this end we shall study the density function 
of a multivariate normal distribution, which is given by 

(4-2) /(x) = — — — i— — — exp ( -ix^A-^x ) , 

^ ^ ■'^ ' (27r)'-/2det(>l) ^ V 2 / ' 

if A is the covariance matrix of the distribution. 

Our first lemma shows that the determinant of any matrix A e M.r{^) is close to 1 

if e is small enough. 

Lemma 4.1. If e < l/(2r) then for any A e Mrie) we have det{A) — 1 + O(e^r^). 
Proof. Let Sr be the set of all permutations a of {1, . . . , r}. Then we have 

(4.3) det(A) = ^ Sgn((7)ai<,(i) • • • ara{r) = 1 + XI Sg"('^)^l<^(l) ' ' ' "'Mr): 

C^Sr C^Sr 
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where 1 denotes the identity permutation. For < A; < r let Sr{k) be the set of 
permutations a & Sr such that the equation a{j) = j has exactly r — k solutions in 
{1, . . . , r}. Then Sr{0) = {1}, Sr{l) = and more generally one has 

\Sr{k)\ < (^^^_^{k- 1)! < r^ for 2 < A; < r. 

Moreover, note that |aio.(i) • • -0^(7(^)1 < e'^, for all a G Sr{k). 
Hence, we deduce 

r r 

sgn((7)ai(,(i) • • • ara{r) = XI 5Z ^g"('^)^i^(i) ■ ■ ■ "-Mr) < ^{(^r)'' < eV^. 

aeSr k=2 aeSrik) k=2 

Inserting this estimate in (4.3) implies the result. □ 

In order to understand the behavior of the density function /(x) we need to deter- 
mine the size of the entries {djk} of A"-*^, if A e M.r{e). The next lemma shows that if 
e is small then the diagonal entries are close to 1 and the off-diagonal ones are small. 

Lemma 4.2. If e < l/(2r) then for any A G Mrie) we have 

l + 0(eV) ifj = k, 
0(e) if J ^ k. 



Proof Recall that 



where Mj^j is the minor of the entry akj which is given by M^j = det{Ai;j) and A^j is 
the matrix obtained from A by deleting the k-th row and the j-th column. 

First, we determine the size of the diagonal entries ajj. In this case, remark that 
Ajj e M.r-i{^)- Hence, it follows from Lemma 4.1 that 



Now, we handle the off-diagonal entries. For 1 < j ^ k < r, let Bj^k denote the set of 
all bijections a from {1, . . . , r} \ {j} to {1, . . . , r} \ {k}. Then, we have 

\Mjk\ ^ \det{Ajk)\ < X Yl \^na{n)\- 

For < Z < r — 1 we define Bj^kQ) to be the set of bijections a e Bj^k such that the 
equation (T{m) — m has exactly r — 1 — I solutions. Since a{k) ^ k then it follows that 
Bj,fe(0) = 0, and more generally one has 

\Bj,k{l)\< (^/_~'^_^il-'^V-<r'-\ forl</<r-l. 
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Hence we obtain 

7 — 1 r— 1 

\Mjk\ < XI n \^na[n)\ < ^r^~^^^ < ^■ 

1=1 aeBj^kH) ^<n+j<r 1=1 

Combining this bound with Lemma 4.1 yield the desired bound ajk <S e. □ 

We know that the Fourier transform of a multivariate Gaussian of covariance matrix 
A is (up to normalization) a multivariate Gaussian of covariance A^^. The last ingre- 
dient we need to prove Theorems 1.1-1.3 is an approximate version of this statement 
when A e Mr{e). 

Lemma 4.3. Let r >2 be a positive integer, R > be a real number and x e M^. 

If e < l/(2r) then for any A e A4r{^) we have 

^ ' Jm<R V 2 ; (27r)'-/Met(A) 2 J 

+ o(e.p(-^!)). 



Proof Since exp (— ^t^At) is the Fourier transform of the multivariate normal distri- 
bution whose density equals 



/(x) 



(27r)'-/2 det(A) 

then the Fourier inversion formula yields 
(4.4) 



exp 



x^^-^x 1 , 



{2ny 



AtlXl+-+trXr) / ^T^^ \ 



(27r)^/2det(>l) 



exp ( — x^A ^x ) . 



Moreover, since \ajk\ < l/(2r) for j k then 



^ ^ ^jktjtk 
l<j¥'k<r 



j=l 



by the Cauchy-Schwarz inequality. This implies 

r r 1 

(4.5) t^^t ^J2Y1 ^iktjh >7^Y. 

Hence, we get 



3=1 k=l 



3=1 



{2nr [ 



exp ( —t^At 1 dt 



>R 



< c^^r [ 
j\\t\ 



>R 



exp — 7||t||^ dt <^ exp ( — 



which in view of (4.4) completes the proof. 



□ 
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Proof of Theorem 1.1. To lighten the notation we shall write 5q for 5q-^ai,...,ar ^-^d Hq for 
IJ-q;ai,...,ar- Let R — 3>/Var(g) logq. First, using Lemma 2.3 we derive 
(4.6) ' 

Sq= dfiqivi, ...,yr)= L^y^^...^y dfiq{yi, ...,yr) + (exp (-2 log^ q)) . 

Jyi>y2> ->yr J \y\^<R ^ 

Next, we apply the Fourier inversion formula to the measure /ig to get 

\yU<R \yU<R •'^^^ 

Since the Fourier transform fiq{si, . . . , s^) is rapidly decreasing, we shall deduce that 
the main contribution to the integral over W of e'*^*^^^"^""'"*''^'''/iq(si, . . . , s^) comes from 
a small ball centered at 0. Indeed, we infer from Proposition 3.2 that 

+ 0(exp (-21og'g)), 

where e = 3(Var(g))~^/^ log g. Hence we get 
(4.7) 

6q = {2nr / e^^^^y^+■■■+^^y^^flq{s^, Sr)dsdy + O (exp (- log' q)) , 

since <C exp (r log g). Now, we make the change of variables 

tj :— •\/Var(g)sj and Xj , for all 1 < j < r 

VVar(g) 

to obtain 
(4.8) 

Sq ^{27t)-^ [ [ ^^it,.,+■■■+trXr)|^ ( ^1 _ _ tr \ ^^^^ 

+ O (exp (-log^g)) . 

Replacing fig I , , . . . , , ) by the approximation (4.1) that we derived in 

\ y/^a.T{q) \/ya.T{q) J 

Proposition 3.1 yields 

^ = i^^r / e^(*-i+-+*-'^) exp i-b^Ct) dtd^ + E,, 

where 

E,<^q-'/'{logqf'- <^q-'/\ 
since (i(g) = g"^^-' and t^Ct > by (4.5). Furthermore, applying Lemma 4.3 we derive 

(4.9) 5g = , , i, [ exp (--^^C-'x) + (q-'/^) . 

^ ' ^ (27rY/^ det(C) h>^2>->xr ^\ 2 J ^ 
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Since Cjk = -Bq(aj, afe)/Var(g) <C (logg)~^ for j 7^ /c by Lemma 2.2, then there exists 
an absolute constant ao > such that C G M.r{l3) with /3 = ao/logq. Therefore, 
appeahng to Lemma 4.2 we obtain 




which foUows from the Cauchy-Schwarz inequahty. Hence we deduce 



(4.10) - 77 1 + llxlr < --x^C-'x < -- 1 - |lx 

^ ^ 2 V loggy " " - 2 - 2 V log?; " 

for some absolute constant cti > 0. This imphes 

L>..>->x ^''P ^"^^ I exp ( -^||x|n dx < exp (-log'g) . 

Inserting this estimate in (4.9) and using Lemma 4.1 we get 

Let Khe a. real number such that \k\ < air/ logg. Since the function ||x||^ is symmetric 

in the variables {xj}i<j<r we obtain 

(4.12) 

exp (-i(l + .)||x|p) 1 (-1= £exp [-{l + ^)y') dy)' 



(27r)-/ 



Xl>X2>--->Xr 



1 f^f 
= - exp O 



r!(l + «;)''/^ r! \ \logg'^ 
The theorem follows upon combining this estimate with (4.10) and (4.11). □ 

Proof of Theorem 1.2. The result can be obtained by proceeding along the same hues 
as the proof of Theorem 1.1, except that we make a different choice of parameters in 
this case. Indeed, choosing it! = 5-\/Var(g')r logr and using Lemma 2.3 and Proposition 
3.2, we obtain analogously to (4.8) 
(4.13) 



+ O (exp (— 4r logr)) . 
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Moreover, we infer from (4.1) that 



(4.14) 5, = (27r)-'- / . ^ . ...^^^ / e^(*-i+ -+*'--'^) exp f -^t^Ct) dtdx + E„ 

i„i / „i„„1. l|t||<31ogq V / 



where 
(4.15) 

Note that 

1 /■ (lOVrlogr)'" / rlogr ' 
ax = — / dx = ; = exp h Ofrloglogr) , 

Xl>X2> ->Xr /l I n t' V 2 V & & V / ' 

by Stirhng's formula. On the other hand, it follows from (4.5) that 

1 -r. \ , If ( W\ , 1 



^ exp ( -^t^Ct ) dt < I exp ( ) dt 



Therefore, inserting these estimates in (4.15) and using the classical bound d{q) — 
exp (0(log 5/ log log q)) we deduce 

/I / logg \\ 
E2 <^ exp — (log q-\- r log r) + O — h r log log r ) ) + exp (— 4r log r) . 

V 2 yiogiogg yy 

Continuing along the same line as in the proof of Theorem 1.1 we obtain analogously 
to (4.11) 



(4.16) .,= (l + o(jll-))^£__exp(-i 



^x^C-^x)cix + £;3, 



where 



<^ exp ( — -(logg + rlogr) + O \ - — h r loglogr ) ) + exp (— 4rlogr) 

V 2 \loglog? 

Furthermore, it follows from (4.10) and (4.12) that 

exp I — -x^C~^x I dx = — exp I O 



(27r)^/2 A,>.,>...>.. ^ V 2 ; r! V Vlog? 



3;i>a;2>--->a:r 

2 



exp I — r log r + r + O ( log r + 



logg 



by Stirling's formula. Inserting this estimate in (4.16) completes the proof. 

□ 

Proof of Theorem 1.3. Since iXq-^ai,...,ar is absolutely continuous with respect to the Lebesgue 
measure, it follows from (1.1) that 

^q;a\,...,a, — i ^q;ar,a\,...,a, — i ~l~ ^q;ai,ar,--;a, — i ~l~ ' ' ' ~l~ ^q;ai,...,a, — 1,0^- 
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Hence, if 2 < s < r < (f){q) are positive integers then 

(4.17) max <5g;ai,...,a. < max 5g.bi,...,bs- 

(ai,...,ar)&Ar{q) {bi,...,bs)£As{q) 

On the other hand, using Theorem 1.2 with s = [(1 — e/2) logg/ log log g], we get 
max Sg.i,i,...,bs = exp ( -s log s + s + O flog s + ^ -j^. 

{bi,...,bs)€As{g) V V logg// g^"' 

The theorem follows upon combining this inequality with (4.17). □ 
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